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ON SOME POLITICAL AND LEGAL QUESTIONS RELATIVE TO OUTER 
SPACE 


MARION MUSZKAT 
Faculties of Law and Social Science, Hebrew University, Tel-Aviv Branch 


The view that war as an instrument of world politics is becoming outdated because 
of the danger to the very existence of the entire human race with which the employ- 
ment of modern armies is fraught, would appear no longer to be questioned. 


Yet there are enormous political, economic, social, psychological and other difficul- 
ties that stand in the way of implementing the idea of disarmament as a means of 
obviating the factors most threatening to world peace. One of these is undoubtedly 
connected with the answer to the question recently posed by Seymour Harris in the 
“New York Times”, viz.- “Can we prosper without arms?” 


Most of the relevant findings indicate an affirmative answer to this question. The 
tremendous possibilities open to advanced nations, particularly the great powers 
(who are, incidentally, the largest producers of armaments), to assist underdeveloped 
countries to find their place in modern civilization, are pointed out. Moreover, it 
is considered no more than the bounden duty of such nations to do so, while it is 
noted that the satisfying of the needs involved, on a commercial basis, constitutes 
a gigantic, many-sided undertaking and, what is more, 2 not unprofitable one. 


If one takes into consideration the fact that the probable cost of construction of 
a single manned space vehicle designed to land on the moon, for instance, and return 
to the earth’s surface, is likely to be rather more than 20 billion dollars, one can say 
that the numeros categories of scientists and technicians concerned in research and 
industry connected with the exploration of outer space may well be proud of the fact 
that their work points to the existence of new possibilities of utilizing, profitably, 
the means at present absorbed by military budgets. This fact probably gives rise to 
one of the most important political and economic aspects of the manifold problems 
connected with outer space. 

These problems without doubt contribute to the growing conviction that a) inter- 
national cooperation, especially scientific and economic, ought to be continually 
extended if mankind is to be spared a great deal of unnecessary effort; b) the reduction 
of the elements which hinder such cooperation and which also constitute a danger to 
peace, can and should be gradually effected; and c) the recent agreement, for instance, 
prohibiting the use of Antartica for any military purposes whatsoever, can and must 
be followed by an agreement outlawing the use of outer space for military ends and 
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establishing international machinery for supervising the adherence to such an agree- 
ment. 

It may thus be clear why, for various reasons, outer space has become the concern 
of the United Nations as well. 

Para. I(d) or the Resolution of the General Assembly of the 13th December, 1958, 
adopted at the Assembly’s 792nd plenary meeting gave rise to the appointment of an 
ad hoc committee for the Peaceful Uses of Outer Space, which decided that in order 
to fulfil its mandate it has 1) to select and define problems that have arisen or are 
likely to arise in the near future, in the carrying out of space programmes, 2) to divide 
the problems into two groups, those which may be emanable to early treatment and 
those which do not yet appear to be ripe for solution, and 3) to indicate, without 
definite recommendation, various means by which answers to such problems might 
be pursued. 

It should be pointed out that the U.S.S.R. and certain other states decided to remain 
aloof from the work of this Committee, as they were of the opinion that the establish- 
ment of the Committee discriminated against the states of the eastern bloc. But it is 
apparent that it was not the struggle for enforcing the principle of equal representation 
of both political camps in the Committee which was the only or main reason for the 
U.S.S. R.’s absence from this body. 

It should be borne in mind that the United Nations resolution establishing the 
Committee arose out of two varying approaches to the question: one, based on the 
assumption that attention should be paid to the use of outer space for peaceful and 
scientific ends only as conceived by the United States, had already found expression 
in the American proposal laid before the General Assembly on the 14th January 1957, 
and the other, supported by the U.S. S. R., pointed out the necessity for considering 
the problem of excluding outer space from the province of military activities as part 
of the general disarmament plan. At the end of 1959, however, during the 14th Ses- 
sion of the United Nations General Assembly, no doubt under the influence of the 
factors favouring the already months old relaxation of international tension, a reso- 
lution, jointly supported by all states, was adopted establishing a new United Nations 
outer space Commission. 

This Commission will be composed of 24 members, 12 representatives of the 
Western bloc, 7 of the Eastern bloc, and 5 representatives of states not in either of 
the two political camps. It was also resolved to call a widely based international con- 
ference to discuss all problems connected with the peaceful uses of outer space. As 
evidence of the growing universal cooperation in this sphere, we may also mention the 
agreement reached in December 1959 at The Hague by the Committee on Space 
Research (COSPAR) established in October 1958 on the initiative of the International 
Council of Scientific Unions (I. C. S. U. ) with a view to ensuring the cooperation of 
all scientific institutions dealing with research of relevant problems, and to convening 


international conferences to discuss questions concerning the exploration of outer 
space. 
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As regards the identification and listing of legal problems which might arise in the 
carrying out of programmes to explore and exploit outer space, the work done by the 
ad hoc Committee will most probably serve as a basis for the relevant discussions in 
the new United Nations Committee on this matter. 

From the results of the work referred to above the following are the legal problems 
susceptible of priority treatment: (1) The question of the freedom of outer space for 
exploration and use; (2) Liability for injury or damage caused by space vehicles; 
(3) Allocation of radio frequencies; (4) Avoidance of interference between space 
vehicles, aircraft etc. 

1. As regards the first question, the Committee took into account the premise of 
the permissibility of launching space vehicles and their flight, once launched, without 
protest from any state, regardless of the territory “‘over’’ which they pass during their 
flight through outer space, and noted the factual recognition of a rule to the effect 
that outer space is freely and equally available for exploration and use by all in ac- 
cordance with existing practice, and that this rule might find its expression in custo- 
mary law or in appropriate future agreements. 

Like the freedom of the high seas, this rule may be termed the freedom of outer 
space, or freedom of the universe, its exploration and exploitation. But it would be 
of paramount importance to bring about an agreement limiting this freedom to the 
peaceful uses of outer space and outlawing its exlpoitation for military ends. It is 
thus regrettable that the Committee did not find it necessary to formulate any proposal 
on this subject, if one bears in mind that its terms of reference were to ensure the use 
of space for exclusively peaceful purposes. 

2. With respect to the second question, the Committee was of the opinion that 
early consideration should be given to the conclusion of an agreement for the compul- 
sory submission of disputes between states as to their liability for injury or damage 
caused by space vehicles, to the jurisdiction of the International Court of Justice. 
The Committee also pointed out that as far as liability for surface damage caused by 
space vehicles is concerned, the Rome Convention of 1952 on Damage Caused by 
Foreign Aircraft to Third Parties on the Surface and the International Civil Aviation 
Organisation’s experience could be taken into account. 

The problem no doubt involves the elaboration of international standards regarding 
safety and precautionary measures governing the launching and control of space 
vehicles. There is also some room for doubt if the experience of the I. C. A. O. as to 
the kinds of injury for which recovery may be had is entirely appropriate. The ques- 
tion, for instance, is whether there should be liability, without regard to fault, for 
certain or all activities, or whether liability should be based on fault. Again, should 
different principles apply depending on whether the injury occurs on the earth’s sur- 
face, in the air space or in outer space; should liability arising out of launching be 
unlimited in amount, or joint or several, where more than one state or an international 
agency participate in any particular activity. It is just because of the differences be- 
tween aircraft and space vehicles that the kind and extend of damage caused by the 
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latter cannot be compared with the effects of damage and injury resulting from the 
use of aircraft. 

A possible approach to the solution of the above problems lies not only in studying 
analogous material based on existing arrangements, but also in coming to original 
conclusions as a result of a joint effort on the part of lawyers as well as experts, 
properly to evaluate the new phenomena and their possible consequences. Difficulties 
may arise, as a consequence of the Committee’s proposition that disputes be com- 
pulsorily submitted to the jurisdiction of the International Court of Justice, as some 
states may prefer other methods of settling disputes. 

3. As to the third question, it was rightly emphasized that to avoid paralyzing 
interference with radio transmissions, a rational allocation of frequencies for com- 
munication with and between space vehicles is imperative, and that the relevant re- 
commendations of the international Telecommunications Union and its special 
committees should be taken into consideration. 

4. As regards the fourth question, technical studies by interested governments with, 
if necessary, the assistance of competent specialized agencies, may have to be under- 
taken in order to prepare the basis for agreements with respect to the prevention of 
physical interference between space vehicles, particularly rockets, and conventional 
aircraft. 

Questions of identification and registration of space vehicles and of the co-ordi- 
nation of launchings, the return to earth and landing of space vehicles are also among 
the problems concerning outer space which are of primary significance. 

Identification of space vehicles can be effected by reaching agreement on the allo- 
cation of individual call signs, orbital or transit characteristics and suitable markings. 
Such an agreement will become increasingly necessary and urgent as the number of 
space vehicles increases and particularly when they begin to return to earth. 

There can be no doubt that identification would be facilitated by a system of re- 
gistration of the launchings of space vehicles, their call signs and current orbital and 
transit characteristics. Registration would help to avoid the potential overloading of 
tracking facilities, and enable other states to distinguish between the various space 
vehicles registered and other objects and to take appropriate measures, if necessary, 
to protect their interests. 

Identification, registration and co-ordination of launchings also constitute the most 
important elements in the solution of problems concerning property, traffic security, 
rescue proceedings, liability and jurisdiction. 

It would seem that these problems could, in large measure, be solved by analogy 
with the rules governing the identification, registration and co-ordination of the 
various kinds of craft, and particularly with those governing navigation co-ordination 
and especially flights of conventional aircraft. 

As regards artificial earth satelites, however, and the other types of rockets launched 
for exploration purposes only, it ought to be recognised that there should be no 
pretensions to property in these objects which should be considered as “res nullius” 
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-rather like a ship abandoned by its crew, a wreck, or a sea-bottle. Unmanned space 
vehicles and such vehicles when manned in the future can not yet be considered as 
“legibus solutus” in any case. It would be most appropriate, for various reasons, to 
bestow national status on space vehicles. This would, first of all, ensure order and 
resolve questions of liability for activities in space, and for the return to earth and 
landing of space vehicles. It may some-times become necessary for space vehicles 
returning to earth and landing to pass through foreign air space or land on foreign 
territory, whether this is occasioned by accident, mistake or distress. Such even- 
tualities make necessary the conclusion of a multilateral agreement governing, inter 
alia, the means of returning space vehicles to the state which has launched them and, 
in the future, their personnel as well, when they are manned. 

The United Nations ad hoc Outer Space Committee was of the opinion that the 
rules of International Law already existing with respect to aircraft and airmen landing 
on foreign territory through inadvertence, accident or distress, might be applied to 
this field as well. 

According to the Committee, certain problems, such as the question of determining 
where outer space begins, in relation to the question of protection of public health 
and safety; the question of safeguarding against contamination of outer space or 
contamination from outer space, in relation to the exploration of celestial bodies; 
and the avoidance of interference of space vehicles with one another, do not appear 
to be ripe for solution. 


As regards the first of the above problems, it may be affirmed that in accordance 
with existing agreements, a state’s sovereign rights come to an end where the air space 
ends. The Roman rule “‘cuius est solum, eius est usque ad coelum”’ is, of course, no 
longer applicable. According to this rule, the moon might be considered as under 
the jurisdiction of the state in the morning and subject to that of another in the 
evening. There is, however, no consensus of opinion among experts as to the location 
of the limits of sovereignty of states in the space above the land and water over 
which they exercise sovereignty, and the proposition that the determination of such 
limits be based on the physical characteristics of the atmosphere and of aircraft, is 
fraught with great difficulty. Demarcation on the basis of present experience and 
knowledge may soon become out of date. The U. N. ad hoc Committee therefore 
came to the conclusion that this question be left to concretization by custom which 
will prevent the acceptance of a boundary so low as to interfere with existing 
aviation regimes or so high as unreasonably to fetter activities connected with the 
use and exploration of outer space. 

It should, however, be borne in mind that territorial water limits became established 
on the basis of security needs. The 3-mile limit was established hundreds of years ago 
in conformity with the then-existing artillery ranges, and this limit is still preserved 
although today deprived of any military significance. 

Taking current scientific, technical and military knowledge and experience into 
account, it may not be unreasonable mechanically to establish by multilateral agree- 


6 MARION MUSZKAT Bull. Res. Counc. of Israel 


ment the altitude limits of the exclusivity of the rights of states. This might be done 
in order to avoid interference which might hinder scientific exploration above the 
earth’s surface and at the same time to regulate the use of space for military purposes, 
for the freedom of outer space, unlike the freedom of the high seas, should and can 
be limited to the right of all states of peaceful exploration and exploitation. 

It is self-evident that the demarcation of outer space for the purpose of outlawing 
its use for military ends would constitute a valuable contribution to the advance of 
disarmament programmes. 

If such an approach could be accepted as suitable for prompt implementation with 
the concomitant setting up of a supervisory body, speculation as to the need of any- 
thing similar to territorial waters or such as a contiguous zone on the border between 
the atmosphere and outer space, would lose all significance. 


With regard to questions relating to the exploration of celestial bodies, the con- 
ventional methods of acquiring territory by declaring new-found areas to be under 
the sovereignty of the particular state which has discovered them etc., would appear 
to be outmoded, for the simple reason that the challenge of the universe can most 
probably be met, to a large degree, by common effort of a real international nature 
undertaken by international organizations in order to ensure to all nations the even- 
tual benefits of an enterprise, being the common achievement and the expression of 
the common will of mankind, which will come to fruition with the advance of science 
and technology. 


The acceptance of a rule that celestial! bodies are incapable of appropriation to 
national sovereignty and that their exploration and exploitation be carried out ex- 
clusively for the benefit of mankind as a whole and under international administration 
ought to be followed by the establishment of another rule to the effect that any acti- 
vities of the above nature should be so conducted as not to destroy any life pheno- 
mena on such bodies. 

It is clear from what has been said above that many complex political and legal 
problems which have arisen as a result of the advance of astronautical science and 
technology demand the almost urgent creation of a new branch of law to govern them. 

Whatever this law may be called - astronautical law, the law of the universe, 
interplanetary law, or the law of outer space - it already exists in an embryonic form 
and has now to be progressively but rapidly developed in order to keep up with the 
needs of reality and the requirements relating to the strengthening of peaceful co- 
operation and prevention of war. 

Based in a certain measure on analogy with the provisions of maritime and air law, 
but at the same time having to provide original solutions for new problems, the new 
law must be ready to deal in the not so distant future with questions which are bound 
to arise when space travel is no longer a mere challenge to future generations. 


COSMIC RADIATION PROBLEMS IN ASTRONAUTICS 


Kurt SITTE 
Department of Physics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


Recent progress in the interpretation of intensity variations of cosmic radiation has de- 
monstrated the existence of streams of high-energy particles emitted during solar activity, 
and of magnetic clouds within the solar system which trap and scatter cosmic ray particles. 
For further study of these phenomena, research is being carried out in the Physics Depart- 
ment along three main lines. By operating meson detectors at sea level and underground, 
data on the energy distribution of the solar component are collected. By observing D-layer 
reflection with an ionospheric recorder—which according to theoretical considerations 
should be strongly affected by cosmic radiation—in correlation with solar activity, the 
low-energy end of the solar emission is checked. In addition, correlation of intensity 
variations with the observations of the solar physics laboratory is expected to yield further 
information about the nature of solar activity. Preliminary results of these experiments 
will be presented. 


I. As soon as space travel became a distinct possibility for the near future, and hence 
the object of serious research instead of a mere dream, it also became obvious that 
the problems of cosmic rays in relation to astronautics, their possible effects on a 
space ship and its crew, had to be studied seriously. After all, we are dealing with 
an enormously penetrating ionizing radiation, by far exceeding in range and in 
complexity of its interactions any known terrestrial source, and we are well aware 
of the dangers of exposure even to the milder radiations of our laboratories to men 
and materials. Complete, or even partially efficient shielding is impossible; it was, 
therefore, imperative to know precisely the amount of radiation and radiation 
damage we have to reckon with when we venture to leave for an extended period 
the protection offered by the earth’s atmosphere. A priori, it is quite conceivable 
that it might prove a prohibitive dose and hence would make manned space flight 
an impossibility unless or until we discover completely new ways to overcome this 
obstacle. 

Fortunately, early investigations gave very encouraging results. When our first 
vehicles, balloons and rocket sonds, reached what might be called the rim of the 
atmosphere—that is to say, atmospheric depths of a few g/cm? ora fraction of a per- 
cent of the entire air layer— the total intensity of the cosmic radiation was found 
to be much less than what we consider in our laboratories a “permissible dose”’. 
However, two new developments have shaken this optimism in recent years, and 
have brought it home to us that much more work will have to be done before 
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we can say with confidence that man is safe from radiation danger in space. I refer 
to the observations which confirmed the existences of vast regions of magnetic 
clouds in the galaxy, and to the discovery of the Van Allen belts’ of radiation en- 
circling the earth. 

As to the first, it is easy to see that such clouds of ionized matter containing, 
and carrying along, a magnetic field will act as traps for high-energy particles. The 
density of radiation inside a cloud will, consequently, differ from that outside: in 
other words, the cosmic ray intensity is by no means quite uniform throughout the 
galaxy or even large sections of it, but may well show a “patchiness” of considerabie 
degree. If that is so, how can we be certain that the intensity measured in the im- 
mediate vicinity of the earth is a true measure of the total radiation exposure ex- 
pected in space flights, or even a fair average? It is not possible that we live in a 
sort of enclosure of low radiation density, but may encounter much higher doses 
once we leave its safe boundaries? 

Perhaps such a conjecture seems unduly pessimistic. It is, however, not without 
some foundation in facts. Or better, there are some facts known to us which may 
well yield exactly this result, although we do not yet have a clear enough picture 
to give a conclusive evaluation. Recent observations of solar phenomena have 
demonstrated that the sun is continuously emitting streams of particles, possibly 
all along its equatorial region. This emission would, as it were, carry away a solar 
magnetic field (if the sun, like the earth, has a more or less permanent, dipole-like 
field). We might view it as a regular wind in that tenuous atmosphere of the sun— 
to which we shall have to refer still often—blowing radially with an intensity fluc- 
tuating with the solar 11-years’ period, and at times erupting into short-lived hurri- 
canes: the solar flares. In this picture, the observed slight periodical variations of 
cosmic ray intensities find their explanation, and likewise the occasional sudden 
intensity fluctuations which are as a rule accompanied by other disturbances like 
magnetic storms. But unless we assume that at least during part of these periods 
the solar wind does not reach us, that the effective range of the sun’s activity ends 
within the earth’s orbit, it also compels us to admit the serious imperfection of 
our knowledge of the primary cosmic ray intensity at some distance from the earth. 

The discovery of the Van Allen radiation belts was even more of a shock, and 
for a while really seemed to make radiation damage a serious threat to space flight 
instead of a remotely possible danger. Perhaps one should have expected it; yet 
everybody was surprised when in the early satellite flights it was found that the 
cosmic ray intensity—supposed to be constant once we have left the earth’s atmo- 
sphere behind—rather suddenly began to increase to a manifold of its supposed 
“saturation value”, reaching amounts which would make long-term continuous 
exposure impossible. However, it was soon learned that this belt does not extend 
very far, and that beyond it the cosmic ray intensity drops back to its normal value. 
Or almost: for further study revealed the existence of a second Van Allen belt, 
still very much more densely filled with corpuscular radiation and very much wider 
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in its extension both in latitude and in distance. The two belts do in fact present a 
danger which is by no means negligible, and it is easy to predict that in the near 
future, a great deal of work will be done in order to elucidate origin and nature of 
the Van Allen belts. At present, we believe we have a consistent theory of the second, 
outer, radiation belt, but the views on the inner belt are still conflicting. A little 
more about this will be said later. 


II. Turning now to the very modest contribution the Technion Physics Depart- 
ment can make to these efforts, I wish to describe three approaches from which we 
hope to obtain results of some significance. From all-too obvious reasons, we cannot 
indulge in the luxury of sending equipment to satellite or lunik heights; and while 
we may perhaps entertain reasonable expectation of flying some apparatus to moder- 
ate heights (as rocket heights go), the bulk of our equipment is earth-bound. Indeed, 
part of one of the experiments of astronautical interest is carried out underground, 
in a room just off the tunnel of the Carmelit underground railway. Let me begin 
with discussing this experiment. 


In this investigation, w-meson telescopes are used to extend the study of time 
variations of cosmic radiation over a large energy region. The connection with 
astrophysical problems follows from the general description of the current views 
onthe causes of these variations given above. Without goinginto details, it will be 
clear that the various theories hitherto worked out (e.g. Morrison”, Parker?, Elliot* 
either predict precisely the energy dependence of the phenomena and of the degree 
of their correlation, or can be tested by observation of such a dependence where 
none, or no regularity, follows from the theoretical assumptions. This is, after all 
nothing more than a conclusion from the obvious fact that any magnetic trapping 
mechanism can operate efficiently only up to energies E < 300 Hp, where H stands 
for the average field intensity and p for the dimension of the trapping field. 


In all the work hitherto reported, the energy dependence of events like flares or 
slow intensity decreases (“‘Forbush effect’’) has only been studied from a comparative 
survey of observations at different stations, i.e. on the basis of the geomagnetical 
latitude effect. This puts an upper limit of about 15 GeV to the primary energy in- 
volved, since particles of higher energies have too high a magnetic rigidity to be 
affected by the earth’s field. It would, however, be of importance to learn about the 
presence or absence of a solar component at energies above 15 GeV, and about the 
possibility of trapping such particles. In order to do this, we have constructed tele- 
scopes of different opening angle, and have put one such telescope underground. 


(Incidentally, the use of the word “‘telescope”, long accepted in cosmic ray work, 
will seem appropriate also to others, now that astronomy with invisible radiation 
has become so important a part of this science. We are now beginning to learn to 
decode all kinds of messages sent from the stars, not only those emitted in visible 
light, and radio or cosmic ray telescopes take their place along with the “classical” 
optical instrument). 
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The principle of the experiment is demonstrated in Figure 1 which shows a schem- 
atic drawing of the two basic units in use. Each consists of four trays of counters, 
and is subdivided so that both threefold coincidences of the lower trays, and four- 
fold coincidences of the entire system are continuously recorded. The particular 
geometry chosen for the telescopes makes the two telescopes ABC and A’B’C’D’ 
“cubical” (i.e. the distances AC and A’D’ equal the linear dimensions if the square 
trays), while the threefold telescope A’B’C’ is ‘“‘semi-cubical”’ (i.e. the distance A’C’ 
is equal to one-half of the tray dimensions). ABCD forms a “‘narrow-angle telescope”’, 
admitting only particles of near-vertical incidence. The dashed lines in the figure 
indicate the average zenith angle of the particles registered in the three telescopes: 
obviously the largest for the semi-cubical arrangement. It will be noticed that by 
recording simultaneously all three rates, we do in fact carry out observations over 
five distinct average inclinations: By taking the difference between the counting 
rates of, say, the narrow-angle and the cubical telescope, we record only those 
particles which pass ABC without traversing D, and hence come with a well-de- 
fined zenith angle between the two others. 

Now the average energy of the mw-mesons arriving with different inclinations 
depends in a well-known way on this zenith angle @; in first approximation it is 
proportional to cos~!0. Wecan, therefore, ascribe to all our five rates corresponding 
average energies: the lowest to the narrow-angle telescope, the largest to the semi- 


Figure 1 


Schematic diagram of the counter arrangement: Narrow-an : i 
lagram : -angle telescope ABCD; cubical telescopes 
ABC and A’B’C’D’; semi-cubical telescope A’B’C’—The dashed lines indicate the average ment 
angles of the three arrangements. 
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cubical setup. But having defined the average meson energy, we can also determine 
the average energy of the primaries which originated these particles. To use again 
for simplicity first-approximation expressions, we notice that in the energy region 
involved (~50-100 GeV) the number of secondaries n, in a nuclear collision in- 
creases with the primary energy Ey according to n,« E)'/* so that we may use the 
following relation between the average energies E, recorded, and the desired primary 
energies: 


B= kz” 


Consider now the underground telescope. It is of the semi-cubical-cubical variety, 
defining the same zenith angles as its prototype above ground but demanding pene- 
tration of an additonal 30 m.w.e. Writing ho for the effective ‘‘absorber thickness”’ 
of the atmosphere, A for the amount of material traversed in addition to the 
atmosphere, and 0; for the average zenith angle of the i-th telescope, we get for the 
average primary energy E%) of the particles originating the registered u-mesons 

EO = k ho th Ns 
ho cos 0; 


h = o for the stations above ground). Figure 2 shows the relative positions, in the 


N(>Eo) in arbitrary units 


: 10! 10° 
Eo x10" ev 
Figure 2 
Primary cosmic ray spectrum showing the average primary energies of the mesons in the five groups 


recorded above ground (open circles), and in the three groups of the underground measurements 
(full circles). 
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well-known cosmic ray primary spectrum, of the five average energies defined by 
the three telescopes above ground, and of the three groups recorded underground. 
A dashed line indicates the upper limit of the region which can be studied by geo- 
magnetic comparison: it is evident that much new ground is covered in this ex- 
periment. 

Unfortunately, cosmic ray data can be collected only slowly, and it would be 
most unwise to jump to conclusions from a still insufficient body of information. 
Just as an illustration, I may mention that after well over a year of continuous 
operation, we still do not feel that we have enough material to compute with satis- 
factory accuracy our own set of coefficients for reduction of the atmospheric effects: 
our statistics, we feel, is still not good enough. We are, therefore, not yet ready to 
give final results on other and more interesting phenomena either. However, two 
examples may be presented which will demonstrate qualitatively the value of the 
experiment. In Figure 3 and Figure 4 I show the records of the narrow-angle, the 
cubical, and the semi-cubical telescopes during two periods of some 10-14 days. 
The first was a period of somewhat higher than normal registered cosmic ray in- 
tensity but of no particular solar activity. We are, therefore, not certain whether 
the observed variation may be purely of atmospheric origin, or may be due to primary 
fluctuations. The second period represents the “tail end” of a prolonged Forbush 
decrease, and hence is certainly due to solar activity. In both cases six-hourly aver- 
ages are plotted in logarithmical scale, so that the relative variation can be compared 
immediately. We note that in Figure 3 no serious and systematic difference is found 
between the three curves, thus indicating terrestrial origin of the variation. In Figure 4, 
on the other hand, the differences are rather conspicuous and demonstrate by the 
smaller effect on the the semi-cubical telescope the decreasing influence of the 
variation-producing agency on primaries of higher energy. The experiment is being 
continued, but it will take at least another year before we shall be ready to give 
detailed results. 

A second study, likewise only in its beginning stage, attempts to add information 
on low-energy primaries to the observations on the high-energy region discussed 
above. It is based on hitherto unpublished calculations by Segal® on the contri- 
bution to the ionisation of the upper atmosphere due to cosmic radiation. His re- 
sults are shownin Figure 5 together with data of some other authors referring to the 
contribution of other factors. It is seen that in the D-layer region, cosmic rays are 
a very important, perhaps the decisive, contributor. But this effect will, of course, 
impose no demands on particle energies—other than those following from the geo- 
magnetic cutoff. Its study—the study of variations and fluctuations in the D-layer 
ionization—is therefore an extension to a few GeV of the work already mentioned. 
However, in order to disentangle the various components, it will also be necessary 
to observe secular variations on ‘“‘quiet” periods, and to compare results with those 
of stations on other geomagnetic location. Some advance in this direction has re- 
cently been reported by Bailey®, but much more remains to be done. 
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Figure 3 


6-hourly rates of the narrow-angle, the cubical, and the semi-cubical telescope during the period 
May 15-27, 1959, in logarithmical scale. 
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Figure 4 


6-hourly rates of the narrow angle, the cubical, and the semi-cubical telescope during the period 
August 6-18, 1959, in logarithmical scale. 
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Figure 5 


Number of ions per cubic centimeter produced by a number of typical sources; (1) Ly, — line, 
(2) Lg — line, (3) L — continuum, (4) total electromagnetic solar radiation, x —solar x-rays, 
(5) cosmic radiation. 


The equipment used in our experiment is basically a model C-4 ionospheric re- 
corder, kindly loaned to us by the National Bureau of Standards in Washington. 
As it is well known, this apparatus emits a powerful pulse of short radio waves, 
rapidly sweeping through a large frequency region. By timing the echo signal, the 
distance of layers of a given electron density is then determined. An example of a 
typical echo sweep is shown in Figure 6. 


Unfortunately, we have not been able recently to operate the recorder at full 
power, since for its new location the antenna has to be built. This is a job which we 
cannot handle with our own technical staff, and we have not yet succeeded in getting 
the responsible authorities to act. However, we hope that this will be possible in the 
very near future, so that we shall be in a position to correlate ionospheric—i.e. low- 
energy cosmic ray—data with the recordings of meson telescopes. 


With the ionospheric we are reaching what might be called the outskirts of the 
Van Allen belt, and we thus are reminded again of one of the most important out- 
standing questions: that of the origin of the inner belt. Two radically different 
views have been offered: according to one (Singer) it results from the decay of 
albedo neutrons, and thus is relatively stable and time-invariant. According to the 
other, however (Morrison®), the inner belt is not basically different from the outer, 
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but is like it the result of the capture by the earth of streams of particles emitted 
by the sun. If this is so, it becomes likely that the inner belt is not constant in time— 
just as the outer is not—and may even vanish at times, undergoing in general variations 
and fluctuations in correlation with solar activity. 


It is clear from this brief description that simultaneous observations of cosmic 
ray intensity variations, solar activity, and ionospheric conditions (not only in the 


Figure 6 
Typical photograph received with the ionospheric recorder, Model C-4. 


D-layer, but especially in the highest layers close to the Van Allen belt) are of con- 
siderable interest. For this purpose, we plan to put to use an optical telescope equipped 
with an alpha filter in addition to the two types of apparatus already described. 
Visual observations will give us information on the state of solar activity, and in 
some cases a warning before an event, and will be supplemented by the other methods. 
Unfortunately, here again no systematic results can be reported as yet, since only 
now the observatory building is at last in construction. Up to this time, we could 
operate the telescope only by looking out of a laboratory window, and hence for 
short intervals which made regular observations impossible. 


Thus, I am afraid that this report contains a few ideas, a lot of promises, and 
scarcely any results—none if we ask for definite conclusions. But we must take the 
data as they come—we cannot manufacture them; and we must work under the 
conditions we meet — we cannot ignore them. This may or may not be an excuse: 
but it is an apology, the only one I have to offer for the paucity of this report—and 
the last one I intend to offer before presenting a fully documented account of our 


work, perhaps a year from now. 


CISION SS 


KURT SITTE Bull. Res. Council of Israel 


REFERENCES 


VAN ALLEN, McILWAIN AND Lupwia, 1959, J. Geophys. Res., 64, 271; VAN ALLEN, J. A. AND 
FRANK, L. A., 1959, Nature, 183, 430. 

Morrison, P., 1956, Phys. Rev., 101, 1397. 

Parker, E. N., 1956, Phys. Rev., 103, 1518; Meyer, P., PARKER, E. N. AND Simpson, J. A., 
1958, Phys. Rev. 104, 768. 

Etuiort, H., private communication (in print). 

SEGAL, I., unpublished private communication. 

BAILey, D. K., 1959, Proc. Inst. Rad. Eng., 47, 255. 

SINGER, S. F., 1959, Phys. Rev. Letters, 1, 171; 3, 188. 

Morrison, P., private communication (in print). 


A PRACTICAL METHOD FOR CALCULATION OF PRESSURE 
DISTRIBUTION ON AIRFOILS IN SUPERSONIC ROTATIONAL FLOW* 


A. KOGAN 


Department of Aeronautical Engineering, Technion-Israel Institute of Technology, 


Haifa 


ABSTRACT 


NOMENCLATURE 


eu =w7) 


= pressure coefficient 


airfoil profile 
undisturbed flow Mach number 
pressure 

gas constant 

entropy 

temperature 

cartesian coordinates 
velocity components in 
x and y directions, in 
units of limiting 

speed 


En 
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Vo,00 = 


A transformation of coordinates, which introduces Crocco’s stream function as independent 
variable, is used to develop a method for calculation of pressure distribution at the surface 
of an airfoil in a supersonic rotational flow. The preparation of tables for the practical 
application of the method is discussed in detail. 
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shock wave profile in (€,7)- 
plane 


= slope of shock wave 
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(Ky,/Kw)o = ratio of shock wave curvature to airfoil curvature ar leading edge. 
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w = condition at airfoil surface; — = condition of undisturbed flow; , = condition 
behind shock wave at leading edge 
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A method for the evaluation of a second approximation of the pressure coefficient 
at the surface of an airfoil in a supersonic rotational flow was developed in reference 1. 
Its main features will be briefly outlined here for ease of reference. The method is 
based on an application of Crocco’s srteam function. The basic equations of flow 
were transformed by the introduction of the abscissa € =x and Crocco’s stream 
function 1 = (x,y) asindependent variables. The ordinate y = y(é,n) satisfies the 
following differential equation: 


yiQ-1) (y+1)/(y- 1) 


u(c? = u?)yeg— (9 = Duv(L—w) yy + “uw? - w’) Yai 
1 wW@-1) 
= ates =s6*)(l-— w7) ce 
as dn (1) 
Its derivatives are connected with the velocity components by the equations 
[ aa 
| alae? 
Piiga eeane 
q u (1 2. a ii (2) 
The tangency condition at the airfoil surface is given by 
y(6,0) = g (6) (3) 


where g(¢) is the known airfoil profile. The boundary conditions at the shock wave 
surface are expressed by 


(4) 


The values of u,v y, and y, behind the shock front can be calcuiated by Eqs. (2) 
and (4) and tabulated as functions of the shock wave slope. 

If 1 = ‘¥(€) represents the shock wave profile in the (.7) plane then d'¥/d. can 
also be calculated and tabulated as a function of 7: 


ia (5) 
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The functions y(€, 7), u(é, n), o(€ ) and S(n) are developed in power series of n 


YolS) + yi(E)n + yo(E)n? + «.. 
uo(E) + u,(E)n + u(y? +4 .. 
v = vo(8) + 0,(6)n + v2(S)n 7+ ... 
S == 9S. Sin + -Soy” +22 
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II 


(6) 


Then the pressure coefficient at the airfoil surface is given by 
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Here p,/p, and p,/p, are the pressure and density ratios across the shock wave 
at the leading edge. The value wo, of the velocity at the airfoil surface is obtained 
from its value at the shock front by use of expansions (6) 


2 2. 2 
( we = u2 + 02 


(8) 


Uo =u — uy = Uyn” a neers 


2 
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Thus in order to obtain a first approximation to c, we neglect all terms O(7) in 
the above expansions. The velocity components at the airfoil surface are approxim- 
ated by their values behind the shock wave for the same value of ¢: 


Up =U; Vg =v (9) 


These are known as functions of the shock wave parameter t. In order to establish 
their functional dependence upon € it is observed that to this approximation 


yon) © Yo'(S) = 9'(S) (10) 


For a given airfoil geometry and angle of attack g(€) is known. The tabulated values 
of u, v and y, can be used in connection with Eq. (10) in order to obtain c,(¢) to a 


first approximation. 
A second approximation to c,(€) is calculated by retaining the first order terms 


in y in Eqs. (8) 
Up =U— UN; Vo =V— 4H (11) 
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The first order coefficients u, and v, are obtained from the expansion of Eqs. (2) 


by oN u{2c*y,(é) a UDYnVi(S)} + 0(7) 
Vile tes) (12) 
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Here the coefficient y,(é) is determined from the first approximation to the 
differential equation (1): 
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The expressions yo"(€) and y,'(€) appearing in Eqs. (12) and (13) can be obtained 
by numerical differentiation of the tabulated values of y.= yo (€)+0(y) and 
Vn = ¥x(6) + O(n), which are related to € by Eq. (10). Similarly 4 in Eqs. (11) is ob- 
tained by numerical integration of Eq. (5), the tabulated values of t, y- and y, being 
related to € through Eq. (10). 

The parameter S,/c, in Eq. (13) is a constant for each undisturbed flow and 
airfoil geometry and orientation. It is calculated from known flow parameters 
behind the shock wave at the leading edge: 


Si she 4y(y cae 1) [(M{ ne 1) Ge Ty 1]? g"(0) (Ky,/Ky)o (14) 
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Finally the dependence of t and of the other tabulated parameters on € is iterated 
so as to satisfy again the flow tangency condition at the airfoil surface: 


v— 04n 
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= 9 (8) (16) 
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In reference 1 a number of numerical examples, covering a wide range of rotational 
supersonic flows past sharp airfoils, have been worked out numerically and the 
results were compared with those obtained by the method of characteristics. The 
present method proved quite reliable and accurate. One of its inherent advantages 
consists in the fact that most of the numerical work involved in its application 
could be performed in a general way, independent of the particular airfoil to which 
it will be applied. Flow parameters like u, v, Yes VY» etc., behind the shock front 
can be tabulated once and for all as functions of one common shock wave parameter, 
e.g. t. The availability of such tables would greatly simplify the numerical work 
involved in the application of the method to a particular flow configuration. 


The preparation of such numerical tables has been undertaken. They will be cal- 
culated by the WEIZAC electronic computer of the Weizman Institute of Science. 
We shall here discuss some of the details concerning their construction. 


Eq. (7) may be written in abbreviated form 
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T) being the value of t at the leading edge. 
The expression wo in Eq. (17) can be developed by substitution from Eqs. (11) 
to (15) into the following form: 
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The following 10 parameters must be calculated and tabulated: 


A, B,C. D 
E 

Er 

u 

Ves Vn 

Av /dé 


(Eqs. 21) 
(Eqs: 22523) 
(Eq. 18) 
(Eq. 4) 

(Eq. 2) 

(Eq. 5) 


At first sight it would seem natural to choose t as the independent variable and to 
calculate these 10 parameters as functions of t at constant intervals At. Such a scheme 
would not be convenient, however, since it would have complicated the procedure 
to be followed in the use of the Tables. As seen from Eqs. (20) the calculation of 
cy along a given airfoil involves numerical differentiation of y’(é) and y,(é and 


the numerical integration 
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Should the Tables be computed at constant intervals of t, these differentiation 
and integration processes would become quite cumbersome since they would have 
to proceed numerically at varying intervals of €. 

Yet there is, of course, no possibility of constructing the tables at constant inter- 
vals of ¢ since the functional relation t = 1(€) is determined by the profile and angle 
of attack of each individual airfoil. 

The practical way out of this difficulty is to calculate the Tables at constant inter- 
vals of yz. The differentiations and integration will then be performed with constant 
increments of yz according to the formulae: 


FAG ie 
dé/dyo 
" dy, 
V6) 379 (= — yi 
1 dyo e) 
yo (5) 
7 = i oD dyo 


% dE ¥O 


The application of the tables to the calculation of the second approximation of 
the pressure coefficient at the airfoil surface will require the following steps: 

1. Place origin of (x,y)-plane at airfoil leading edge and draw x-axis in direction 
of undisturbed flow. Apply Eq. (10) to the given airfoil profile g(x) and calculate 
€ =x at constant intervals Ayz of yz. 

2. Apply Eqs. (24) to obtain by numerical differentiation and by numerical inte- 
gration yo "(é), y,'(é) and n(€), respectively, at constant intervals Ay, of yz. 

3. Use the Tables to obtain u, A, B, C and D at constant intervals Aye. 

4. Calculate g"(0); read up from the Tables the value of E at the airfoil leading 
edge — E(to). 

5. With the values obtained by steps 2, 3 and 4 calculate uo, v9 and Wg at constant 
intervals Ayz, by Eqs. (19) and (20). 

6. Calculate (1 — wy7)*” 

7. Obtain from the Tables the value of F at the airfoil leading edge — F(to). 

8. Calculate by Eq. (17) c, at constant intervals Ay,. 

9. Insert values of uo and vg, as obtained by step (5), into Eq. (16), to obtain 
the improved correspondence between the calculated values of c,, step 8, and the 


abscissae x = €. 
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ABSTRACT 


A study is made of the flow past a supersonic airfoil which oscillates harmonically at 
high frequency. An asymptotic solution of the linearized equations of motion is ob- 
tained by a method analogous to the “‘eikonal’’ representation of optical waves. The 
boundary conditions are satisfied on the surface of the airfoil, so that the effects of airfoil 
shape are accounted for. The results include (1) the first approximation (for high frequencies) 
to the velocity potential, (2) the second approximation to the oscillatory pressure on airfoil 
surface. 


EQUATIONS OF MOTION 


The supersonic flow past a harmonically oscillating airfoil is governed, within the 
linear approximation, by the partial differential equation’ 


a*p* a7p* 


+ 210 2e" = 2g =0 (1) 


2 — 
Fie 2) Ox? oy? Ox 


In eq. (1), as in the following, #*e'' is the oscillatory velocity potential, w denot- 
ing the angular frequency and t the time; x and y are dimensionless Cartesian co- 
ordinates, stationary with respect to the airfoil and measured in units of the airfoil 
chord c, the x axis being parallel to the undisturbed flow; M denotes the free-stream 
Mach number and 4 is the dimensionless ‘‘reduced’’ frequency defined by 


MC 
ao 


A= 


where dy denotes the speed of sound. The boundary condition which ¢* has to 
satisfy on the airfoil surface can be stated in the general form 


(ee) goes cW (2) 


on n=0 


Received January 17, 1960 
25 Bull. Res, Counc. of Israel, Vol. 8C, 1960 


26 MEIR HANIN Bull. Res. Counc. of Israel 


where n denotes the (dimensionless) normal distance from the surface of the airfoil, 
and We’ is the normal component of the oscillatory velocity of the air on the 
surface. When the shape and the motion of the airfoil are given, the values of W on 
its surface are known. In addition to (1) and (2) the velocity potential has to satisfy 
the “‘radiation’’ condition! according to which ¢*e'®’ must behave at large dis- 
tances from the airfoil like an outgoing wave of finite energy. 

In the following, an asymptotic representation of the velocity potential will be 
obtained for large values of the reduced frequency 4. For this purpose we put 


OM = FO Veer (3) 


and we assume that F does not vary rapidly with position as 4 tends to infinity. 
Thus G will determine primarily the phase, and F the amplitude, of the oscillatory 
flow; the surfaces G(x,y) = const will represent the wave fronts. An analogous 
assumption is made in theoretical optics* to account for the transition from wave 
theory to geometrical optics; the function G is called there the ‘‘eikonal’’. The 
assumption permits us to expand F in an asymptotic series of powers of 47: 


F= AF (a tFOG,y) + 4? FO (yy) +... 1. (4) 
Substituting (3) and (4) into the differential equation (1) and separating out the 


resulting powers of 4~* we find that G, F“, F‘ etc. must satisfy the partial dif- 
ferential equations 


ee (M? —1) (5) +2M—— =1, (5) 

By ay * (MOP 0] Fa 5 (FG a9 2] = 0,6 
ay ay [Oe] es gr 

it Fae ECV IAT) he (for m>2) (7) 


Similarly, substituting (3) and (4) into the boundary condition (2) we find that 


k=0 (8) 
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and that G, F, F™etc. must satisfy on airfoil surface the conditions 


(G),=0 = 9, (9) 
CO iy ee (10) 

ana ego 
(ele ean (ie Ayres i 
on he 1 an ra or m2). (11) 


The equations (5)-(11) show that the required asymptotic expansion of the 
velocity potential can be obtained by the following procedure. First the phase function 
G(x,y) is evaluated by solving the equation (5) under the condition (9). When G 
is known, we evaluate the function F‘” by solving equation (6) under the condition 
(10), and so on. It should be noted that the ‘‘radiation’’ condition is not taken 
into account in the determination of the solution; however it will be shown that 
the resulting solution satisfies the ‘‘radiation’’ condition automatically. 

As the boundary conditions (9)-(11) are prescribed on the airfoil surface, the 
solution will depend on the shape of the airfoil. We represent the surface of the 
airfoil by the parametric equations 


X=X(0), Y= Yo(9) (12) 


it being assumed that the parameter o increases along the surface in clockwise 
direction. We also define the (dimensionless) length /(¢) along the surface by 


1'(0) = [(xo'(0))? + (oo)? T° (13) 


where the prime denotes differentiation with respect to o. The angle 9(c) between 
the tangent to the surface and the positive x axis will the be given by 


sin 9(¢) = yo'(a)/I'(a), cos 9(a) = X9'(a)/I'(a) (14) 


the tangent being assumed to point in the direction of increasing a. The local (di- 
mensionless) radius of curvature R(a) of the airfoil surface, 


R(o) = —1'(0)/3'@) (15) 


will also appear in the results. 


THE PHASE FUNCTION 


We turn now to the evaluation of the phase function G(x,y). Applying to equation 
(5) the usual method of characteristics for first-order partial differential equa- 


tions’ we denote 
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0G 0G 


bets 2 2 
ee oe saat ree H = — +2Mp—(M*-1 | (16) 


so that the differential equation (5) becomes 


H=0: (17) 


According to the method of characetristics, the solution of (17) is obtained by 
solving the system of ordinary differential equations 


dx dy dG dp 
Ae TEP hy oer iad Ripa i te PH coin a) 
d 
<= (Ho + H,) (18) 


where the subscripts indicate partial differentiation. In our case this system reduces to 


dx dy dG 


FW se 2u% gee hdiide Jaks a 2525 2 2 
eae (M* — 1)p, aan Mp (M1) pg, 

dp dq 

‘ds : oe” 1) 


The quantity s introduced in (18) or (19) is a parameter which varies along the char- 
acteristic lines (“‘rays’’) of the solution. The initial conditions for the system (19) 
are obtained by taking 


s=0 on airfoil surface; (20) 
the boundary condition (9) gives then 
G=0 when s=0 (Bid i) 
and from (12) we have 


X=xX,(0), y=yolo) when s=0. (22) 


The values po(o), qo(o) which p and q respectively assume on the airfoil surface 
must satisfy 


H (p= poo), 4 = qo(a)) = 0, (23) 
Po(F) Xo (F) + Gol) Yo'(a) = 9 (24) 


(the ‘‘strip’’ relation (24) expresses the fact that G does not vary along the surface). 


The equations (23) and (24) determine po(o) and qo(c); with the aid of (16) and (14) 
we find after some calculation that 
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sin 9(c) cos X(c) 
Opies See es ae < 
Polo) hes Vache ae I — Msin Xo) a) 
The conditions (21) and (22), together with the statement 
P= Pol), ¢=4() when s=0 (26) 


provide a complete set of initial conditions for the differential equations (19). Inte- 
grating (19) under these initial conditions we obtain the solution 


rh M — sin Xo) > cos (a) 
X =X(o) + 1=M sino) Ss, y=Yo(o) ay — M sin (a) : 


(27) 


1 
~ 1 — Msin Xo)” 


G (28) 
It should be pointed out here that the airfoil under consideration must have a sharp 
leading edge and be sufficiently thin to satisfy 


1 — Msin &(c) > 0 (29) 


along its surface; otherwise a detached shock wave would be formed and the lin- 
earized equations of motion would not be valid. 


The results (27) and (28) constitute the required solution for the phase function 
G(x,y). In fact, from eqs. (27) relations of the form s = s(x,y), o = o(x,y) can in 
principle be obtained, which, when substituted into (28), would give the function 
G(x,y). Alternately, G(x,y) can be found by eliminating o from 


X =X9(a) + [M — sin 8(o)]G, y = yo(a) + cos X(a)G. (30) 


The nature of the solution is, however, understood best when (s,o) are regarded 
as independent ‘‘characteristic’’ coordinates which define the points (x,y) through 
the equations (27). When a is kept constant and s varies from zero to infinity, the 
equations (27) represent a straight line (‘‘ray’’) emerging from the point (x9(¢), Vo(e)) 
of the airfoil surface. Because of the restriction (29), the ray lies entirely in the field 
of flow (so that it does not cross the region occupied by the airfoil). As o varies, 
the equations (27) define a system of rays, each ray emerging from a (different) point 
of the surface. Thus the coordinate o identifies the various rays while the coordinate s 
determines the position of a point on its ray. The solution (28) gives then the values 
of the phase function G on the rays. We notice that G grows linearly from zero 
to infinity as the distance along a ray increases. 

The wave fronts G(x,y) = const and the direction of the rays can be found by a 
simple graphical construction (Figure 1). We select a value for G, say G,; then 
from a point O on the surface we draw the normal OQ with length G, and from 
the end Q of the normal we draw the length MG, in the direction of the x-axis, thus 
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Figure 1 


reaching the point P. It follows now from the solution (27), (28) that the point P 
lies on the ray which emerges from O, and that the value of the function G at the 
point P is equal to G,. Furthermore, if we draw PS perpendicular to OQ (i.e. parallel 
to the surface at O) so that the point S lies on the normal OQ, then the length OS 
is equal to the characteristic coordinate s corresponding to the point P. In fact, 
denoting by (xo, Yo) and (x,y) the Cartesian Coordinates of O and P respectively, 
we have from the construction in Figure 1 


Y—Yo=G,cos9, x—x9 =MG, —G,sin9 


and OS = G, — MG,sin9, so that G, = a ee OS 
A comparison of these relations with the solution (27), (28) shows at once that 
the point (x,y) lies on the ray emerging from (X9g, yo) and that 


OS =s, G(x,y)= G,. 


The wave front G(x, y) = G, can therefore be found by carrying out the construc- | 
tion described above for several points (x9, yo) of the surface and by drawing a | 
curve through the resulting points P. It can be proved also that the wave front | 
G(x,y) = G, is the envelope of circles of radius G, with centers at the points 
(Xo(o) + MG,, yo(o)). 

The foregoing solution for the phase function has a simple physical interpretation: 
it shows that the propagation of disturbances from a surface element of the airfoil 
is formed by vectorial superposition of two motions, of which one is the propagation 
with the speed of sound dy directed normally to the surface element, while the other 
is the undisturbed flow moving downstream with the speed Mag. Thus the disturb- 
ance waves are carried downstream by the main supersonic flow while they propagate 
away from the airfoil with the speed of sound. It should be noted that the character 
of the phase function agrees fully with the requirements of the ‘“‘radiation’”’ condition. 
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VELOCITY POTENTIAL 


Having found the phase function G, we proceed to evaluate the amplitude function 
F™ by solving the differential equation (6) under the boundary condition (10). 
For this purpose it is most convenient to employ the characteristic coordinates 
(s,o) as defined in (27) (or in Figure 1). Differentiating the equations and inverting 
the resulting relations we obtain 


oo cos 3 ‘ Cue ee 

see an ee eS 

pera | aes EN ; 

oo M — sin9 Rae cs 

es ee. I SI 

- 7) [1 — Msin9 + qa ; (31) 
se i : M-—sin3 s aes 
ae — | sint? (1 — Msin 3) + Tan R(o) [1 — Msin 3 + RG 4 
eS ean [1 - Msin9 hls Rb ian BINT A 
dy a i—Msin 9 RG) Bee ariel Co 


where 8 = 9c), R(a) and /'(c) are defined in (13)-(15). The relations (31), (32) allow 
us to express in terms of s and o the partial derivatives which appear in the dif- 
ferential equation and the boundary condition for F‘’’. With the aid of the solution 
(28) for G, we get 

CRD FORM 


6G OF? aG 


= Py» oO = ; 
éy = oy t [M asl? a Ox ae By 
0*G E a°G 1 
dy? ee ioe (34) 

[1 — M sin 9(o)] R(o) +s 
0G OG 0G 1 
a SF ae ae Sa 35 
on Ox Sika oy oo83 1 — Msin 9) 22) 
Therefore the equation (6) and the condition (10) become 

ory 1 1 (1) 
= ae Le at a f= 0, 36) 
Os i 2 [1 — M sin 9(¢)] Ro) + ( 
(F) 9) = —icW(o) [1 — Msin 4(0)]; (37) 


in (37) the notation W(a) is introduced to indicate that the normal velocity component 
W on the surface of the airfoil is a known function of position along the surface. 
The equation (36) shows that the partial differential equation for F’» reduces to 
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an ordinary differential equation along the rays; this is true also with respect to the 
other functions F“”” since (33) remains valid when F“? is replaced by F“™. It follows 
that the motion and the geometry of each surface element of the airfoil determine 
completely the high-frequency flow on the entire ray which emerges from the element. 


Integrating (36) under the boundary condition (37) we obtain the required ampli- 
tude function F": 
S rf 
[1 — M sin 9(c)] R(o) 


FY = —icW(o)[1 — Msin 9(0)] [1 + (38) 


This result, together with (4), (8) and (28), gives the first approximation to the velocity 
potential for high frequencies: 


‘ = is 


[1 — M sin 3(a)] ma €@ 1—Msin Go) ihe 
+0(A~*) 


o* = —icW(o)2 [1 — Msin 9@)] [1 + 


The solution (38) shows that the amplitude F“ of the velocity potential decreases 
to zero as the distance from the airfoil along a ray tends to infinity, provided that 
the radius of curvature R(o) is finite at the point from where the ray emerges. When 
a surface element is flat so that its radius of curvature is infinite, the amplitude 
F™ remains constant on the corresponding ray. Therefore, the velocity potential 


given by (39) represents a wave of finite energy, thus satisfying fully the ‘‘radiation”’ 
condition. 


The higher terms F” of the asymptotic expansion (4) can in principle be found 
in the same way as F“"?. Unfortunately the non-homogeneous term in the differential 
equation (7) turns out to be very complicated, so that higher approximations to the 
velocity potential in the flow field will not be considered here. We notice, on the 
other hand, that the boundary condition (11) allows us to find the values of F on 
the airfoil surface directly from the first approximation F“ given above. This enab- 
les us to evaluate the pressure on the surface within second approximation. 


PRESSURE ON AIRFOIL SURFACE 


We denote by P*e'™ the oscillatory pressure on the surface of the airfoil. The linear- 
ized Bernoulli equation for unsteady compressible flow gives 


pt = — Poo cags + me) (40) 


i" 
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where fo is the density of the undisturbed air. Substituting for the velocity potential 
p* the high-frequency asymptotic expansion 


g* a Bite ate, a 72 FY ‘a ov *)] e 4G (41) 


(eqs. (3), (4) and (8)) and recalling that on the airfoil surface G = 0 and s = 0, we 
find that the pressure amplitude P* on the surface is given by 


SE an ere Gln) oF”? : 
pt = — Pore (i Ma )FO 42 [ia - MSF +M=—| +00, ”)| 


Ox s=o (42) 


All the terms appearing in this equation can be evaluated by employing the results 
obtained above. From the solution (28) and the relations (31), (32) it follows that 


0G 1 
a ieee peo oe 
Ox 1 — Msin9 (43) 


while the equations (11) (with m = 2) and (35) give for F on the surface 


OF 
ee eh ge =i = 
(F?), _, = — i(1 —Msin 9) (— Ne 
OF FO 
(lie i io i 44 
i(d Msin 9) ( ay cos 9 a7 Sin 9) (44) 


When (43) and (44) are substituted into (42), and the partial derivatives of F are 
expressed (with the aid of (31) and (32)) in terms of the characteristic coordinates 
we get 


(1) +p) 
eee an O40 ATs Ola )eao -1¢, _ . OF -2 45 
pe = — For lig at 4 (1 Msin 9) ( x ) {Fo )| (45) 


The values of F‘ and dF’/ds at s = 0 are provided by the solution (38). It follows 

that the oscillatory pressure on the airfoil surface is given, within the second approx- 
imation for high frequencies, by the simple result 

1—Msin9,_, ee 

a j —____—— } 46 

P* = — poagW[1 + i——"* 41 +007) (46) 

where the normal velocity component W, the radius of curvature R and the inclination 

angle 9 refer to the same point on the airfoil surface as the pressure P*. The geometric 

parameters of the airfoil influence the pressure on its surface mainly through the 

second term of (46) which may be important when the reduced frequency A is not 
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extremely large. It is to be noted that the airfoil shape affects also the distribution 
of the normal velocity component W on the surface. In the special case of a flat 
airfoil with zero thickness we have R = oo so that the formula (46) reduces then 
to its first term; this term is found to coincide with the high-frequency limit of the 
Possio result* for flat-plate airfoils. 
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A HYPERSONIC RAMJET USING A NORMAL DETONATION WAVE 
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ABSTRACT 


The utility of conventional ramjets at hypersonic speeds is seriously limited due to perform- 
ance deterioration and growing structural complexity, as flight speed increases. 

Some of the associated problems can be alleviated by operating the engine combustion 
chamber at supersonic mean inlet velocities. The governing equations for engine perform- 
ance with a normal detonation wave in the combustion chamber are derived, and the 
flight regime for which combustion through normal detonation is possible is investigated. 

Performance of an ideal detonative ramjet using acetylene fuel is presented for flight 


Mach Numbers of 6 to 10. 


A ramjet engine configuration is suggested for performance with normal detonation. 


NOMENCLATURE 

C = exhaust nozzle velocity i 

coefficient 
F = engine thrust — lbs t = 
fla = fuel-air ratio (p= 
g = acceleration of gravity v = 
H ~~ = total enthalpy Oe 
M = Mach number yn = 
M, = flight Mach Number w = 
Mp, = detonation Mach Number, y = 

defined by eq. 14 y = 
Dp = static pressure Na = 
tg = total pressure p = 
R = gas constant 
Stations 


0 = free stream stagnation 

2 = entry to detonation wave 
3 = exit from detonation wave 
4 =) nozzle. exit 
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specific fuel consumption — 
Ibs/hr/Ib 

static temperature 

total temperature 

velocity 

flight velocity 

detonation velocity 

engine weight flow — lbs/sec 
ratio of specific heats 
average ratio of specific heats 
inlet efficiency 

density 
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INTRODUCTION 


Considerable attention has been given recently to the use of air-breathing engines 
in missions normally reserved for rocket-powered vehicles. Aside from subsonic 
and supersonic cruising missions, in which the air-breathing engine’s superiority is 
unchallenged, the possibility of using them for booster or hypersonic cruise missions 
merits further investigation because of their low fuel consumption. Both of these 
missions require engine operation at hypersonic speeds, a region as yet unexplored 
by air-breathing engines. 

As is well known current material limitations restrict the maximum allowable 
turbine entry temperature for turbojet engines to about 2000°F. Due to this limitation 
the turbojet cycle is less efficient than the ramjet cycle at speeds above Mach 4!. 
Therefore the appropriate air-breathing powerplant for hypersonic speeds is the 
ramjet, which on occasion has been characterized as a “flying stove-pipe’”’. However, 
it is questionable whether the simplicity and low weight generally associated with 
ramjet engines can be maintained by conventional ramjets designed for hypersonic 
speeds. The following are some of the principal difficulties encountered: 


a) Inlet 


For subsonic combustor performance, the ratio of inlet capture area to combustor 
area rises sharply with flight Mach Number. The capture area ratio for a perfect 
inlet is shown in Figure 1. For attainable inlet efficiencies this ratio would, of course, 
be lower; but, as can be readily seen, the inlet geometry requirements become very 
severe at high speeds. If reasonably good inlet performance is required at lower 
speeds also, say during vehicle acceleration, a variable geometry inlet becomes a 
necessity. Due to the very high air temperatures encountered at hypersonic speeds, 
some of the inlet surfaces will require cooling. In light of these considerations, it 
can be seen that the inlet will constitute a heavy and complex engine component. 


b) Exhaust nozzle 
At hypersonic speeds the exhaust nozzle operates at pressure ratios above 100. 
Exit to throat area ratios required for full expansion are shown in Figure 2. The 
severe geometry requirements are compounded by the high gas temperatures in the 


nozzle. As for the inlet, good performance over a range of flight speeds requires 
considerable variation in nozzle geometry. 


c) Combustion chamber 


Since combustor inlet temperatures exceed 1500°F, the chamber requires regenerative 
cooling. In addition, the high chamber temperatures place a limit on the possibility 
of converting chemical fuel energy into thermal energy of the gas stream. This limit 
is set by dissociation and ionization reactions in the gas stream at high temperatures. 
As a result, progressively less of the fuel energy is available for raising the thermal 
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Inlet capture area to combustion chamber area ratio. 


energy of the gas stream. Unless the energy expended on dissociation and ionization 
can be reconverted to thermal energy by reassociation reactions as the gas cools 
in the exhaust nozzle, it makes no contribution to the propulsive thrust of the engine. 
Although the degree of reassociation to be expected in the combustion products of 
fuel-air mixtures is not completely determinate at this time, there seems little doubt 
that high-temperature dissociation and ionization places an upper limit on the speed 
at which useful performance can be obtained from conventional ramjets. 


THE DETONATION RAMJET 


The difficulties associated with conventional ramjet operation at hypersonic speeds 
can be alleviated to some extent if the combustion process is carried out at supersonic 
gas velocity. Little or no diffusion would then be required prior to combustion, the 
exhaust nozzle would not require a convergent section, and gas temperatures through- 
out the engine would be considerably lower. 

Conventional ramjets have been restricted to low subsonic combustor inlet velocities, 
due to the very low velocities of flame propagation associated with the deflagration 
of fuel-air mixtures. However, combustion waves moving with supersonic velocities 
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(detonations) were first observed over fifty years ago, and have been intensively 
studied since. These studies were restricted to detonation waves propagating into 
stationary gases. Recently, however, the possibility of stabilizing detonation waves 
in supersonic gas streams has been investigated, and the application of standing 
detonation waves to ramjet propulsion has been suggested” *. This has been followed 
by experimental work in which detonation waves have been successfully stabilized 
in supersonic gas streams* °, 

The standing detonation wave is analogous to a standing shock wave established 
in a supersonic stream, and it has been suggested that it is actually a shock wave 
followed by a subsonic combustion wave°. Although on the basis of aero-thermo- 
dynamic considerations alone an infinite number of supersonic detonation velocities 
are possible for a given combustible mixture of gases, only waves travelling at sonic 
velocity with respect to the burned gas, Chapman-Jouget detonations, have been 
observed in experiments on initially stationary mixtures’. 

In a manner analogous to shock wave stabilization in supersonic diffusers, either 
normal or oblique detonation waves may be stabilized in the combustion chamber. 
For a particular detonable mixture, the total pressure loss (or entropy gain) of the 
combustion process will be minimized for high detonation wave angles with respect 
to the mean flow direction, since they occur at the lowest Mach Number con- 
sistent with Chapman-Jouget detonation. The optimum wave angle for best engine 
efficiency will, however, depend on the inlet efficiency. An estimate of the upper limit 
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of detonative ramjet performance can be obtained by assuming normal Chap- 


man-Jouget detonation in the combustion chamber and perfect component per- 
formance. 


EQUATIONS FOR NORMAL DETONATION RAMJET 


Considering the station designations given in Figure 3, the engine exhaust velocity 
V4, iS given by 


DETONATION WAVE 


0 eae 3 4 
Figure 3 
Station Designation for detonative ramjet 


d Yas1 Pag7t 
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where 734 Signifies the average specific heat ratio between stations 3 and 4. 
Assuming one-dimensional flow, perfect gas behaviour, and neglecting changes in 
molecular weight, the total enthalpy ratio across the detonation wave is obtained 

from the first law of thermodynamics as 
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From the one-dimensional momentum equation connecting stations 2 and 3. 


pee 1+ 2M} (3) 
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which by use of the continuity equation across the wave may be expressed as 
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Substituting equation (4) into equation (2) gives 
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But for a Chapman-Jouget detonation, the Mach Number of the burned gases, 
M;, is one, and equation (5) therefore, takes the form 


Hiss ib (3) (2="') (1 + 72M3) (6) 
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Considering the flow to be adiabatic between the inlet and the combustion chamber, 
equation (6) expresses the total enthalpy rate between station 3 and 0, i.e. H3/Ho. 
The total pressure ratio across the detonation wave is given by 


3 
tg La \Wark 
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where the average ratio of specific heats, y, is to be evaluated between the static and 

stagnation temperatures at the station in question. 
From the one-dimensional momentum equation across the wave 
Pal MM 3 (8) 


P2 1+ 73M3 


and substituting equation (8) into equation (7), we obtain the total pressure ratio 
across the detonation wave as 
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For a Chapman-Jouget wave this takes the form 
ise 

ie 1 sak 

P; _ 1+ 72M} 2 


P, ¥3 +1 


Y2 

: 10) 
Yo—1 ,,2\V2—1 ( 
(1+ 5 M3) 


Vol. 8C, 1960 M. ARENS 41 


If viscous and mixing losses are neglected, equation (10) gives the total pressure 
ratio for the engine combustion chamber. 


Substituting equation (6) and (10) into equation (1) gives the engine exhaust 
velocity as 


Ve = Cy (2) (1 + y2M,’) (v2 —1) Hy 
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In order that the wave be stabilized, the gas velocity entering the detonation wave 
must equal the detonation velocity of the fuel-air mixture, 


M, = pete (12) 


or 


—1 —{ 
M,= Mp, 2S (1+ 35 M3) (13) 


where the detonation Mach Number, Mp, is defined as 


M pe esl (14) 
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Therefore the combustion chamber Mach Number for stabilization of a normal 
Chapman-Jouget detonation wave is given by 


BS Mp, 
Ve ee Ue (15) 
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and substituting this equation into equation (11) gives the following expression for 
the engine exhaust velocity 
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The specific thrust and specfic fuel consumption can be obtained from the exhaust 
velocity by the well known relations 


F _ (1 + fla) v4 = % 


- ; (17) 
sfc = ae (18) 


DETONATIVE RAMJET PERFORMANCE WITH ACETYLENE FUEL 


The ideal performance of a detonative ramjet using acetylene fuel is shown in Figure 
4. In calculating the performance perfect component efficiency was assumed. The 
thermodynamic properties of the fuel—air mixture were assumed identical with those 
of air in chemical equilibrium as given in reference 8. The detonation velocities of 
acetylene—air mixtures were taken from reference 7. 

Figure 5 is a plot of Mo/sfc, a measure of overall engine efficiency, using the min- 
imum sfc for each Mach Number as given in Figure 4. Since the heating value of 
acetylene is of the same order as conventional jet engine fuels, the performance 
shown in Figure 5 may be compared to the known performance of conventional 
ramjets at supersonic speeds. In making this comparison it should be recalled that 
the operating inlet efficiency of the detonative ramjet should be considerably superior 
to that of the conventional ramjet due to the much smaller amount of diffusion 
required. 

It can be seen that the engine efficiency of the detonative ramjet at hypersonic 
speeds exceeds that of conventional ramjets at supersonic speeds. As far as fuel 
consumption is concerned, it may therefore be anticipated that hypersonic cruising 
using a detonation ramjet will be more economical than cruising at supersonic or 
even subsonic speeds using conventional powerplants. 
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Figure 4 
Specific fuel consumption of ideal detonative ramjet using acetylene 


GAEL tts 
Perfect Component 

= Performance 
o| & 
=| 

uw 

7p) 
> 
O 
ie 
uu 
Oo 
te 
we 
WwW 
uJ 
z 
oO 
a 

35 9 10 


FLIGHT MACH NUMBER Mo 


Figure 5 
Optimum engine efficiency of ideal detonative ramjet using acetylene 
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The detonative ramjet, like the conventional ramjet, is inoperative below a mini- 
mum flight speed limit. This limit is the minimum speed at which stable detonation 
can be established for the leanest detonable fuel—air mixture. In other words, the 
detonation velocity must not exceed the maximum flow velocity attainable in the 
engine. Ideally this limit is therefore given by 


[ork =, (19) 
Os 
or 
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The maximum permissible detonation velocity as a function of flight Mach Number 
at 100,000 feet is shown in Figure 6. The detonation velocities of two acetylene—air 
mixtures are also shown in the figure. 
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Figure 6 
Maximum allowable detonation velocity for standing detonation wave 
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ENGINE CONFIGURATION 


A schematic jayout of a ramjet using a normal detonation wave is shown in Figure 7. 
The area variation of the combustion chamber must be sufficiently large so as to 
allow stabilization of the detonation wave at minimum flight speed and maximum 
fuel-air ratio, as well as at maximum flight speed and minimum fuel-air ratio. It 
is to be expected that the detonation wave, much like a normal shock in a supersonic 
diffuser, will seek its own position of equilibrium. Increasing fuel-air ratio at fixed 
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Figure 7 


Detonative ramjet 


flight speed will force the wave towards the exhaust nozzle, whereas increasing 
flight speed at fixed fuel—air ratio wil! tend to suck the wave in the direction of the 
inlet. 

Should the engine be required to operate at flight speeds sufficiently low so that 
stable detonations are either not feasible or detonative performance is poor, it is 
possible to design the engine to be capable of alternating between detonative com- 
bustion and conventional subsonic combustion. By forming a throat in the exhaust 
nozzle a normal shock will be forced into the inlet and subsonic flow established in 
the combustion chamber, allowing conventional operation. 


CONCLUSIONS 


Some of the difficulties to be encountered by conventional ramjets at hypersonic 
speeds can be alleviated by designing the engine for supersonic mean entry velocities 
into the combustion chamber utilizing a stabilized detonation wave. The efficiency 
of such an engine is competitive with that of conventional ramjets at supersonic 
speeds. It is possible to size the engine geometry so as to permit alternating between 
the conventional and detonative mode of operation. 
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SOME DEVELOPMENTS OF THE END-BURNING CHARGE 


ANSELM YARON 
Scientific Department, Ministry of Defence, Government of Israel 


INTRODUCTION 


The magnitude of thrust obtainable in the usual type of the end-burning charge is 
limited, once the type of propellant and grain diameter is given. An attempt to in- 
crease the thrust of an end-burning charge rocket motor by means of increasing 
the burning surface will fail because of the strongly regressive characteristic of the 
first phase of any such modified surface. The two possible burning surface modifi- 
cations, of the concave or convex type recede in general into a surface equal 
approximately to the circular cross section of the grain. 

The purpose of this paper is to indicate some possilbilities of a stable thrust aug- 
mentation of an end-burning grain by means of a special grain design. 


END BURNING GRAIN CHARACTERISTICS 
The usual end-burning grain consists of a propellant cylinder inhibited against 
burning on one of his bases and on the lateral surface exposing the other unin- 
hibited base to burning. The burning front is generated as an envelope of spheres 
whose radii are the instantaneous burning distances (Figure 1). Applying this rule 
of burning front formation to an increased initial surface of an end-burning grain 


inhibitor 


Figure | 
we reach the following conclusion: 

a) A convex surface, generated by the revolution of a curve around the axis of 
the grain transforms itself in the course of burning in a concave-convex flat surface 
with a cup; this surface recedes ultimately into a circular surface equal to the cross- 
section of the grain (Figure 2a). 

b) A concave surface, generated by the revolution of a curve around the axis 
of the grain with angles of intersection with the axis of the grain, transforms itself 
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Figure 3a Figure 3b 


into a concave spheroidal surface and ultimately into the cross-section of the grain 
(Figure 3a). 

The transformation of the convex surfaces during burning is due to the fact that 
the burning front at the point of contact with the inhibitor (restriction) generates a 
toroidal surface which expands until it reaches the axis of the grain. This surface 
is predominant and eliminates any previous forms (Figure 2). 

The transformation of the concave surfaces is due to the fact that any angle (e.g. 
the angle of the generating curve with the axis) serves as a focus generating toroidal 
or spherical surfaces which dominantly transform the surface into a spherical con- 
cave surface approaching ultimately the cross-section of the grain (Figure 3b). 


STABILISATION OF THE BURNING SURFACE 
We shall consider a positive (convex) cone surface burning towards the base of the 
grain (Figure 4a). The points of contact of the surface will generate an enveloped 
toroidal surface unless prevented by inserting between the conical surface and the 
restrictive layer a cylindrical layer of a propellant having a higher burning rate than 
the main grains. In such a case the points of contact between the two propellants 
will not become foci of generation of the receding toroidal surface because the 
faster burning layer will reveal the surface of the slower burning propellant parallel 
to the inhibitor, before the main burning front reaches the point of contact. Figure 
4b shows the formation of the new burning fronts at intervals f: v,, v. indicating the 
two burning rates. In the same drawing we can see that a new cone is generating, 


having a different half angle than the original one. The angle of the new cone is 
obviously 
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Figure 4a Figure 4b 


Figure 4c 


Choosing the half angle of the slow burning propellant cone equal to the ratio of 
the burning rates of the two propellants, we shall obtain a constant burning surface. 
The ratio between the circular grain cross-section, usually serving as the burning 
surface and the increased conical surface (neglecting the contribution of the thin 
fast burning layer) is given for a positive cone by: 


The enclosed graph (Figure 5) shows the dependance of the ratio between circular 
grain section area and the conical developed area and the burning rates ratio. It 
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Figure 5 
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appears obvious that an inverse disposition of the fast and slow burning propellant 
layers will behave as a normal convex cone end burning charge, the slow burning 
propellant acting as a burning inhibitor. A negative cone with an inserted central 
fast burning cylindrical core will behave in a similar manner stabilizing itself at a 
surface of a truncated cone. Keeping the ratio of the inner fast burning core diameter 
to the outer diameter small, we can consider the generated surface as being one of 
a cone (Figure 4c). 


FURTHER APPLICATIONS OF THE METHOD 


The mounting of the described conical charges in a rocket motor case might present 
difficulties in design or produce a prohibitive tail-off at the end of the burning. 

It seems, therefore, desirable to shorten the reaction zone. This can be achieved 
using a cylindrical layer of a fast burning propellant embedded in the slower burning 
charge (Figure 6a). A further shortening of the reaction zone could be obtained 
using a combination of a central core and several concentric cylinders made of a 
fast burning propellant and embedded in the slow burning grain (Figure 6a, 6b). 


fast 


Figure 6a 


Figure 6b 


The extension of this procedure cannot, of course, proceed ad infinitum, the thick- 
ness of the interposed slow burning layers being the limitation. The amount of 
layers depends on the capacity of the propellant to burn in thin layers without 
breaking up as well as on the technological difficulties arising in the manufacture 
of the grains. 

Another employment of the layer technique consists in the possibility of a limited 
thrust programming. To achieve this scope it would be necessary to use layers and 
cores of fast burning propellant of variable thickness. 
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The calculation of the generated surface bases itself on two assumptions: 

a) the burning front generates along the instantancous radii of curvature; the rate 
of change of the radius of curvature with time being the rate of burning (slow pro- 
pellant) and assumed constant. 

b) the burning front of the fast propellant is assumed always perpendicular to 
the x-axis and assumed constant as well. 


—-- : een 


ad 
Vea 


new burning front 


x 
Figure 7 
Then, setting the following nomenclature: 
y=I(xX); the equation of the curve of contact between the fast and slow 


propellants, the fast propellant assumed located between the curve 
and the x-axis, and 
avons first second and third derivatives of the y = f(x) curve with respect 
tomx; 
12)3/2 
p= Cae ee the radius of curvature of the said curve. 


u“ 


wy = 04; the rate of burning of the slow propellant taken with the appropriate 
dt : ; 
sign according to the character of the curve, 
& V2; the rate of burning of the fast propellant assumed constant as well, 
: and finally — 

ae H, the ratio of the two burning rates. 
Vv, 
we obtain, after derivation of the radius of curvature with respect to x: 

dp 3(1 ae he ye rs Wz (1 4 yey ye 

= D a a (1) 
dx y. 


and introducing the two propellant burning rates: 
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Solving at first with respect to y” and y’ we have after separation of variables: 
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Integrating, we obtain: 
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Separating again the variables, we have: 
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The expression (4) gives the generated curve in items of its slope and x. 
The final integration giving the explicit form of y as function of x is obtained 
by quadratures. The area is obtained later by a simple quadrature of 
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(x,, X, being the two reference coordinates of the curve involved). 

It has to be stressed that the final integration of expression (4) has to be done 

by development into series of the integral curve which yields: 
see Im] 20, Cia oes 
{1 = [in 2c, (« + C)P}* 


(4a) 


The exact analytical treatment of the solution might present some difficulties 
especially if the basic curve of contact between the two propellants will present 
discontinuities or singular points. In such a case, as well as in cases where the pre- 
cision of the solution is of limited importance, it appears advisable to recur to the 
graphical solution using calculated radii of curvature as reference points. 


CONCLUSIONS 


The described method of end-burning charge surface augmentation can offer a 
vast game of applications in sustainer rocket motors. The fabrication difficulties 
of concentrical charges vary with the type of propellant and in case of cast type 
composite propellants should not present excessive complications. 


IONIZATION MEASUREMENTS IN SHOCKED GASES 


Y. MANHEIMER-TIMNAT 
Scientific Department, Ministry of Defence, Government of Israel, Tel Aviv 


The study of electrical properties of thermally ionized gases produced by shock 
waves, was initiated by A. R. Kantrowitz and his collaborators! , first at Cornell, 
later at AVCO. 

The problem was thoroughly investigated for the first time by S. C. Lin®, who 
used three different methods for measuring the electrical conductivity in shock- 
heated argon: 1) the conventional probe method; 2) interaction between the ionized 
gas and a magnetic field — which he calls the ‘“‘magnetic experiment”; 3) rotation of 
the gas flow and consequent interaction with a magnetic field (‘‘vortex experiment”’). 
The shock Mach range studied was 4-16 and results obtained were in good agree- 
ment with Spitzer and Haerm’s’ theory in the shock Mach range 8-10, showing 
small deviations for higher Mach numbers. 

Lamb and Lin® measured the electrical conductivity in air for shock Mach 
numbers 10-18 by the magnetic method and obtained results agreeing with theory. 

Petschek and Byron* studied the approach to equilibrium in shock Mach range 
10-20. They used two different experimental methods, measuring: 1) the intensity 
of the continuum radiation emitted by the gas; 2) the electrostatical diffusion po- 
tential. They succeeded to show that, while the initial stage of the ionization process 
is dominated by impurities, electron-ion collisions are the rate determining process 
after the number of electron has reached 10% of its final (equilibrium) value. 

Niblet and Blackman® investigated ionization in air by a photographic method 
using a glass tube in which shock waves were produced by discharging a bank of 
capacitors. They measured the time required for obtaining equilibrium in the 
shock Mach range I1-17. 

The use of microwave methods for measuring electrical properties of ionized 
gases goes back to the war period. Saenger, Bredt and Goercke? in Germany mea- 
sured the electron density in the exhaust gases of a V-2 rocket by the reflection of 
centimetric waves, while Andrew, Axford and Sugdon'® in England found the 
electrical conductivity of rifle muzzle flames by attenuation measurements in the 
3 cm range. Microwave methods were used in the post-war years in investigations of 
electrical properties of flames!!! and of discharges in gases (reviews given by Loeb”? 
and Goldstein'.) 
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Low and Manheimer!® !° developed a microwave absorption method, by which 
results on the electrical properties and the ionization rate of shocked gases in Mach 
range 7-10 were obtained. 

The principle of the method is to send a microwave beam through the shocked 
gas and measure its attenuation as a function of time. Experiments were performed 
with air, argon and nitrogen using 3 cm waves. For shock Mach numbers 8.2 to 
10.2 in air the electron density calculated from attenuation measurements is in good 
agreement with that predicted from shock-wave theory assuming thermal equi- 
librium. The most important contribution to the electron density comes from the 
nitric oxide, NO, whose ionization potential—9.25 ev—is the lowest of all species 
present. The ionization time, measured from the beginning of the attenuation till 
its maximum (equilibrium) value is reached, is short (7-12 microseconds). The product 
of this time by the initial air pressure plotted as function of shock Mach Number 
extends the results obtained by Niblett and Blackman® to a lower Mach range. This 
points to the following mechanism!’ for NO formations: 


0,420; O+N,NO+N; N+0,<NO+0 


since relaxation times involved are of the right order of magnitude, according to 
experimental evidence'®*°, 

The experiments with nitrogen were performed with a gas containing 1/4% of 
oxygen. This small quantity of oxygen was sufficient in order that the dominating 
factor in the ionization process in the Mach range studied (7.4 to 8.8) should again 
be nitric oxide. The maximum electron density, calculated from the experiments, 
was reached in a somewhat longer time (about 15 microseconds) than in air, probably 
due to the lower oxygen content. Agreement with electron density predicted assuming 
thermal equilibrium was fair. A rough estimate of the dissociative recombination 
coefficients of nitric oxide shows order of magnitude agreement with coefficients 
of typical diatomic gases, quoted by Loeb!’. An interesting effect, probably due to 
reaction between the oxygen impurity and the hydrogen driver gas after the contact 
surface has reached the test section, was observed. This agrees with similar obser- 
vations by Clouston, Gaydon and Glass?!. 

The argon experiments (shock Mach range 7.1-8.5) showed agreement with equi- 
librium ionization at the lower end of the range only. This seems to be due to the 
fact that the ionization time of argon (100 microsec. at M = 10, see ref. 4) is longer 
than test time in the shock tube used. The gas contained 0.2°% impurities, which 
may account for the quick ionization onset and the large spread of values at high 
Mach numbers. Oscillations in the microwave power output before beginning of 
the attenuation were observed in the case of argon. This suggests partial omission 
of 3 cm waves. 

The microwave absorption method allows one to obtain information on the elec- 
tron density and ionization rate in gases heated by medium strength shock waves; 
in particular relatively low shock Mach numbers, where other known methods gave 
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unsatisfactory results. may be studied. The attenuation range most suited for this 
method lies between 0.1 and 10 db; measurements are possible at lower attenuation 
levels too, providing special precautions to stabilize the oscillator’s output are taken. 
The above mentioned attenuation values may cover a considerable range of state 
variable of the gas, since it is easy to vary the following parameters within rather 
wide limits: 

a) the effective path of the microwave through the gas, which determines the magni- 
tude of the attenuation coefficient; 

b) the microwave frequency; 

c) shock strength and initial pressure. 
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ON THE METEOROLOGY OF THE PLANET MARS 


J. NEUMANN 
Department of Meteorology, The Hebrew University of Jerusalem 


Of the constituent gases of the Martian atmosphere the only one which has been 
definitely identified is CO,. The presence of water vapour is next to certain, although 
its amounts are likely to be small. O, is also probably present, but again in small 
amounts. N,j is likely to be the chief constituent. 

Optical methods lead to the estimate that atmospheric pressure at the Martian 
surface is ca. 85 mb (ca. 64 mm HgO° C), i.e. higher than the pressure at which human 
blood boils. Atmospheric density at the surface is ca. 1074 g cm~?. Thus pressure 
and density at the Martian surface are about one order of magnitude smaller than 
at the Earth’s surface. As the Martian g is small (38% of the Earth’s g), pressure 
decreases but slowly with height and at a height of ca. 30 km pressure is the same 
as that at the same height in the atmosphere of Earth. 

The general level of air temperature is about 30°C lower than on Earth. The at- 
mosphere failing to provide an effective greenhouse protection, the daily range of 
air temperature is nearly 100°C. Under the sub-solar point, temperature reaches 
about 30°C, while at the winter pole temperature drops to ca. — 60°C. 

In all probability, the adiabatic temperature lapse rate (g/c,) is ca. 3.7°C/km for 
the lower Martian atmosphere. Despite this small value, there is not much evidence 
in favour of convective cloud formation. Mainly because of the sparsity of cloud 
cover the planetary albedo is low: 0.15 against the Earth’s 0.40. One consequence 
of the paucity of clouds is that we can often see the solid surface of Mars. 

There is some observational evidence that the winds are mainly westerly. The 
thermal wind equation indicates that the variation of wind with height is moderate, 
again, mainly because of the small Martian g. On the whole, the dynamics of the 
Martian atmosphere is governed by the same forces as is the case for the atmosphere 
of Earth. 
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HEAT TRANSFER STABILITY ANALYSIS OF SOLID PROPELLANT 
ROCKET MOTORS 


REUEL SHINNAR AND SGAN ALUF MENACHEM DISHON 
Ministry of Defence and Israel Defence Forces Ordnance Corps, resp. 


This paper presents a theoretical investigation of the causes of unstable burning in 
solid propellant rocket motors. A possible mechanism of unstable burning is deve- 
loped. Small perturbation analysis is applied to a set of differential equations des- 
cribing heat transfer in the solid as well as in the gas film near the burning surface. 
Criteria of instability are derived directly from the equations having only the physical 
properties of the propellant as parameters. The assumption of a time lag, necessi- 
tated in other theories, is avoided. Thus, a clearer physical picture of the instability 
situation is obtained and verification by experimental tests made feasible. It is shown 
why an endothermic reaction at the decomposing surface of the solid propellant 
favours stability. The stabilizing influence of turbulators and radial holes drilled 
into the charge is explained. 
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